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ABSTRACT 

| Recent observations of non-thermal X-rays from supernova remnants have been attributed to synchrotron radiation 
^ i ■ from the loss-steepened tail of a non-thermal distribution of electrons accelerated at the remnant blast wave. In the 
fNl ■ test-particle limit of diffusive shock acceleration, in which the energy in shock-accelerated particles is unimportant, the 
0® \ slope of a shock-accelerated power-law is independent of the diffusion coefficient k, and on how k depends on particle 
^ ' energy. However, the maximum energy to which particles can be accelerated depends on the rate of acceleration, and 
O ■ that does depend on the energy-dependence of the diffusion coefficient. If the time to accelerate an electron from 
c | ! thermal energies to energy E » m e c is t(E), and if k oc E@ , then t(E) oc in parallel shocks, and t oc E 2 ~@ in 
Q-f perpendicular shocks. Most work on shock acceleration has made the plausible assumption that k oc r g (where r g is the 
Q particle gyroradius), so that fi - 1 at relativistic energies, implying a particular (wavelength-independent) spectrum 
£j , of MHD turbulence, where Kolmogorov or Kraichnan spectra might be more physically plausible. I derive the /?- 
^2 dependence of the maximum electron energy resulting from limitations due to radiative (synchrotron and inverse- 
^ ■ Compton) losses and to finite remnant age (or size). I then exhibit calculations of synchrotron X-ray spectra, and 
model images, for supernova remnants as a function of /? and compare to earlier /? = 1 results. Spectra can be 
\ considerably altered for j3 < 1, and images are dramatically different for values of (3 corresponding to Kolmogorov 
cd or Kraichnan spectra of turbulence. The predicted images are quite unlike observed remnants, suggesting that the 
turbulence near SNRs is generated by the high-energy particles themselves. 



MICROPHYSICAL PARAMETERS OF STANDARD SHOCK-ACCELERATION THEORY 

Diffusive (first-order Fermi) shock acceleration is presumed to provide the relativistic-electron distributions de- 
manded by radio and, in a few cases, X-ray observations of supernova remnants (SNRs). In the test-particle limit 
in which the fast particles exert no influence on the shock, the spectrum is independent of the value and energy- 
dependence of the diffusion coefficient k. When the energy in accelerated particles becomes non-negligible, however, 
various nonlinear effects appear. Most directly, the preshock gas is pre-accelerated by a shock precursor whose extent 
depends on the diffusion distance rjj = K/u s hock which does depend on particle energies. One observable effect is 
concave spectral curvature if k(E) increases with E (Ellison and Reynolds 1991; Reynolds and Ellison 1992). For 
both test-particle and nonlinear acceleration, the acceleration rate for relativistic electrons depends on the diffusion 
coefficient. 

Macrophysical parameters, those required for hydrodynamic and thermal-shock modeling of SNRs, include the 
remnant age t, shock velocity w s h (ideally w s h(0)> preshock density no, and compression ratio r. (A piece of important 
thermal-shock microphysics concerns non-Coulomb electron heating at shock, in the absence of which T e <sc T, 
initially). 



For standard shock-acceleration theory, we need in addition the upstream magnetic field B\, the obliquity angle 
#Bn between the shock normal and B\, and the diffusion coefficient upstream and downstream, k\ and ki (in general 
k{E, r,t)). In the standard picture (e.g., Blandford and Eichler 1987), in strong shocks (M A ,M S > 10), once particles 
have energies well above thermal, the primary scattering is from resonant MHD waves (wavenumber k oc particle 
gyrofrequency Q: Icr = Q/c). Quasi-linear theory allows the inference of a diffusion coefficient. Assume a wave 
spectrum I(k) = Ak~ n erg cnT 3 (cm^'f 1 (so the energy density at k ~ kl(k)). Then 



where r](E) is the (energy-dependent) "gyrofactor," A m f p - jjr g , and 77 is the fractional magnetic energy in waves 
resonant with particles of that energy. In the so-called "Bohm limit", normally taken to give the minimum physically 
possible diffusion coefficient, 77 = 1 . 

So 77 oc fe/fe)) -1 oc k n £ l . Since Q. = eBc/E for extreme-relativistic particles, r\ oc E l ~ n and k\\ oc rjr g oc 
E 2 ~ n - E@. Then for Kolmogorov turbulence we have n - 5/3 => B = 1/3, while for a Kraichnan turbulent spectrum, 
n - 3/2 => B - 1/2. Note that for constant gyrofactor (as is often assumed), B - 1 => n - 1, i.e.. a "white-noise" 
spectrum with equal energy per decade. 

This describes diffusion along magnetic-field lines. Perpendicular (cross-field) diffusion is generally assumed to 
result from a cross-field displacement of one gyroradius for every scattering along the field: 



For quasi -perpendicular shocks (#Bn near 90°) and large 77, k «; k\\ and acceleration can be much faster than for parallel 
shocks (Jokipii 1987). Even if 77 = 1 (Bohm limit), there is some obliquity dependence, mainly from field compression 
where #B n ~ 90°. But for large 77, we have faster acceleration where the shock is perpendicular but slower where #B n 
drops well below 90°, giving larger azimuthal brightness variations. Large 77 corresponds to weaker turbulence, so 
that # Bn is well-defined. When 77 is near 1, the obliquity-dependence is smaller anyway. 

Radio observations of SNRs provide few constraints on the microphysical parameters, because the power-law 
part of the spectrum contains relatively little information: the diffusion-coefficent dependence is only second order, 
through nonlinear effects (e.g., a slight spectral hardening with energy; Ellison and Reynolds 1991). But the maxi- 
mum energies to which particles can be accelerated depend on acceleration rates which strongly reflect the diffusion 
coefficient. Now all known shell SNRs are fainter in X-rays (thermal or not) than the extrapolation of their radio 
spectrum (14 Galactic remnants, Reynolds and Keohane 1999; 11 LMC remnants, Hendrick and Reynolds 2001). So 
any non-thermal emission results from the cutting-off tail of the electron distribution, giving direct information on the 
maximum energies. Such synchrotron X-ray emission has now been identified in several Galactic remnants, inferred 
from lineless spectra that cannot be produced by thermal processes (SN 1006: Koyama et al. 1995; G347.5-0.3, Slane 
et al. 1999; and several others). 

These remnants can then be used to study various microphysical quantities. The magnitude of the diffusion 
coefficient is constrained by the maximum electron energies observed, and its upstream value in particular is related 
to possible "halo" emission from upstream electrons diffusing ahead of the shock. The energy-dependence of the 
diffusion coefficient can be constrained by studying remnant morphology and detailed spectral shape, compared with 
model calculations such as those to be summarized below. The obliquity dependence of acceleration affects the 
azimuthal variation of X-rays around the shell. The overall compression ratio can be constrained by the halo visibility 
(and thermal modeling). Finally, observed maximum energies also constrain the shock speed and other hydrodynamic 
parameters. 




(1) 




(2) 



Then 
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SYNCHROTRON X-RAY OBSERVATIONS OF REMNANTS 



SN 1006, the prototype object for synchrotron X-rays, has a spectrum that is well fit by a model in which 
acceleration is limited by electron escape above some E max (Reynolds 1996). The fitted characteristic rolloff frequency 



Fig. 1. Left: SN 1006 Chandra image (adapted from Long et al. 2002). Right: RCW 86 (SW corner), Chandra 
image (adapted from Rho et al. 2002). 

of 7 x 10 17 Hz (where the spectrum has dropped a factor ~ 6 below the lower-frequency extrapolation; Long et 
al. 2002), with the known mean magnetic field (from inverse-Compton TeV emission, (B) ~ 10 /igauss; Tanimori et 
al. 1998; Dyer et al. 2001) implies E max = 53 TeV. Furthermore, the known shock velocity implies that the observed 
electron spectral cutoff energy is too low to result from synchrotron losses in this magnetic field (t\/2(E msx ) > 1200 
yr) or finite age (t acce \ < 800 yr). This leaves only a change in diffusive properties for /Imhd S2x 10 17 cm to explain 
the observed rolloff. 

The detailed cylindrically symmetric model of Reynolds (1996; 1998) fits the spectrum well (Dyer et al. 2001), 
but predicts an external halo with brightness 5 - 10% of postshock peak just ahead of the shock, and an extent of 
many arcmin, the diffusion length of the high-energy electrons. However, no halo is seen in the deep, low-background 
Chandra image: any preshock emission < 1% of the postshock maximum on the bright rims. Plane-shock models 
seen edge on predict a lower brightness, ~3%, but still in conflict with the observations. While efficient acceleration 
of particles can increase the overall compression ratio (e.g., Berezhko, Elshin, et al., 1996), this occurs over the 
long diffusion lengthscale, while the immediate density jump at the small-scale thermal subshock actually decreases, 
making the problem worse. This problem is discussed at length in Long et al. (2002) where it is suggested that the 
post-shock magnetic field may actually increase by more than due to simple fluid compression, perhaps from some 
kind of cosmic-ray-driven instability (Lucek and Bell 2000). 

In general, the halo can be made unobservable by making it thin rather than faint. For this option, the diffusion 
length Ki/w s h must be less than the width of the sharp edge, < 10" for the Chandra observations. This would require 
an unreasonably high upstream magnetic field, at least 20 yuG, especially given SN 1006's height of 570 pc above the 
Galactic plane. 

In the larger, older remnant RCW 86, weak X-ray lines in SW corner imply anomalous abundances if the contin- 
uum is thermal (Vink et al. 1997), while energetics problems (Rho et al. 2002) afflict an explanation of the continuum 
as non-thermal bremsstrahlung (Vink et al. 2002). However, if the harder continuum is dominated by synchrotron 
emission, reasonable abundances can be accommodated. (Borkowski et al. 2001; Rho et al. 2002). Chandra data im- 
ply rolloff frequencies v c ~ (7 - 10) x 10 16 Hz, consistent with loss-limited acceleration (shock velocities are known 
from Ho- emission to be 600 - 900 km s" 1 ; Ghavamian et al. 2002). For radiative-loss-limited acceleration, 
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so we need 77 > 1 and/or r > 4 again, even for a perpendicular shock. 



CONSTRAINTS ON ENERGY DEPENDENCE OF THE DIFFUSION COEFFICIENT 

I have generalized the modeling code of Reynolds (1998) to include arbitrary energy-dependence of the diffusion 
coefficient. The models assume a spherical blast wave encountering uniform upstream material containing a uniform 
magnetic field whose sky-plane projection is vertical and which makes an angle (p with respect to the sky plane. 
The remnant dynamics are Sedov; each post-shock fluid element is endowed with an electron distribution N(E) = 
KE' S exp(-E/E max ), where the maximum energy is in general a function of shock speed, remnant age, diffusion 
coefficient, and shock obliquity angle 9r d . This distribution is then evolved downstream subject to radiative and 
adiabatic losses; electrons are assumed to have short enough post-shock mean free paths that they stay confined to 
their original fluid element and share its expansion. Details are given in Reynolds (1998). 

Here we generalize by allowing an energy dependence of the diffusion coefficent, k <x E@. For the relatively 
steep spectra of plasma waves predicted in simple homogeneous turbulence models, such as Kolmogorov turbulence, 
8 < 1, so that A\\ rises more slowly with energy than r g . At some high energy Eh, then, we will have A\\ = r g . Above 
this energy, we must have an effective n = 1 (constant "gyrofactor.") So we confine our attention to E < Eh and write 



Eh 1 i~g \ Eh 

where Kb = r g (Eh)c/3 = Ehc/3eB, the Bohm-limit value at E = Eh- 

Then for a parallel shock, Ogn = 0°, the acceleration time is given by (generalizing from Reynolds 1998) 



r(r+\) c E h IE 
B(r-\)e u \B x \E h 



r acc = n/ _ ^ -^—[ — 1 (6) 



whereas for a perpendicular shock, 9s n - 90°, we have 

2r 



(2-B)(r-\) e u\B x \E h 



2-/3 
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with r the shock compression ratio. Note that t(90°)/t(0°) oc (EjEhf 1 ^ 1 ^ so that if 8 < 1, acceleration to energies 
E <sc Eh is much faster at perpendicular shocks than parallel. 

The result of these effects for supernova-remnant morphology is dramatic (see Figure 2). For a spherical remnant 
expanding into a uniform magnetic-field B\, around the "equator" where 6b,i ~ 90° the electron spectrum will continue 
to much higher energies than near the "poles" (##„ ~ 0°) if 6 < 1. However, the spectrum will not differ much from 
the 8 - 1 case because the emission is dominated by the bright equatorial "belt". 

CONCLUSIONS AND FUTURE WORK 



Because it samples electron energies where the spectrum is cutting off, X-ray synchrotron emission from SNRs 
can provide significant information on the microphysics of shock acceleration (mainly on the diffusion coefficient), 
to a much greater extent than radio synchrotron emission. In the remnant of SN 1006, the shock is almost certainly 
perpendicular at the bright non-thermal rims. Simple predictions for the spectral break due to radiative losses give 
too high a break, as do predictions based on the finite remnant age (or size). For the spectrum to break where it 
is observed, the diffusive properties of the external medium probably have to change on scales above ~ 0.1 pc, so 
that electrons escape much more freely when their gyroradii are comparable to this scale. The absence of an exterior 
synchrotron "halo" probably requires a cross-shock jump in magnetic-field strength considerably greater than given by 
the fluid compression ratio there, suggesting the possibility of some kind of immediate post-shock turbulent or plasma 
amplification of magnetic field. In the considerably older remnant RCW 86, synchrotron emission is important, 
though not dominant; careful modeling of both thermal and non-thermal emission is required to extract information. 
The rather high observed break frequencies require rapid acceleration: perpendicular shocks with r\ > 1 and/or r > 4. 

The modeling of shock acceleration for diffusion in steep wave spectra in the context of a spherical SNR shock 
has shown that quasi-perpendicular favoritism, already present in constant-gyrofactor models, is much more extreme 
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Fig. 2. Four simulations for different values of /?. The magnetic field is assumed to have a vertical projection on 
the sky plane, and to make an angle 4> = 30° with that plane. 



if p < 1 (wave spectra steeper than k\ <x k' 1 ). While predicted spectra aren't much different from the/? = 1 (constant- 
gyrofactor) case, predicted images are strikingly unlike real remnants for /? < 1 . A likely explanation is that the MHD 
waves causing diffusion are self-generated, with k~ l spectrum. In general, however, it is clear that we need better 
thermal and non-thermal models for X-ray emission in SNRs. Such improved models ought to allow the extraction of 
a great deal of important and useful information about the microphysics of shock acceleration. 
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